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Abstract 

We show how macroscopic manifestations of P (and T) symmetry breaking 
can arise in a simple system subject to Aharonov-Bohm interactions. Specifi- 
cahy, we study the conductivity of a gas of charged particles moving through a 
dilute array of flux tubes. The interaction of the electrons with the flux tubes 
is taken to be of a purely Aharonov-Bohm type. We find that the system 
exhibits a non-zero transverse conductivity, i.e., a spontaneous Hall effect. 
This is in contrast with the fact that the cross sections for both scattering 
and bremsstrahlung (soft photon emission) of a single electron from a flux 
tube are invariant under reflections. We argue that the asymmetry in the 
conductivity coefficients arises from many-body effects. On the other hand, 
the transverse conductivity has the same dependence on universal constants 
that appears in the Quantum Hall Effect, a result that we relate to the validity 
of the Mean Field approximation. 
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I. INTRODUCTION 



More than thirty years after it was first studied, the interaction of fiux tubes and charged 
particles -the Aharonov-Bohm (AB) proble continues to show a tantahzing richness. In 
recent years it has received renewed attention due mainly to its intimate relation with Chern- 
Simons (CS) field theory and the ocurrence of fractional statistics in planar physics. One can 
implement intermediate statistics on particles by attaching to them appropriate 'statistical' 
(i.e., fictitious) electric charge and magnetic fiuxB. Two such particles are then subject to 
AB scattering, and acquire a statistical phase when moving around each other. Varying the 
values of the charge and fiux one can continuously interpolate between bosons and fermions. 
Hence their generic name, anyons. 

On the other hand, CS field theory provides a very convenient way to implement AB 
interactions between particles. The CS field strength vanishes outside isolated singularities; 
thus, there are no classical forces among particles. However, it is expected that a mean field 
approximation, in which the particle is regarded as moving in an average uniform magnetic 
field, should be valid provided the fiux per particle is small and that classical trajectories 
enclose a large number of fiux tubes. This approach has been succesfuUy employed to 
calculate various properties of the anyon gasi'u. In Ref. ^ a self-consistency argument is 
given to support the mean field approximation. 

An important issue in fractional statistics and AB systems is that of discrete (parity 
and time reversal) symmetry breaking. This problem becomes of notorious relevance when 
one considers the possible presence of anyonic excitations in high-Tc superconductors, which 
are believed to present time reversal non-invariance. Having these ideas in mind, March- 
Russell and Wilczekl sought for manifestations of P and T violation in CS theories. They 
found that the scattering cross section of identical anyons (which is directly related to the 
AB cross section) shows an asymmetry provided an additional, non- statistical interaction 
is included. The parity violating term arises from the interference between the amplitudes 
of AB scattering and additional interaction. Caenepeel and MacKenzieH have used these 
results to test further the validity of the mean field approximation. They consider the 
motion of one anyon as the incoherent sum of two-particle scattering processes. They find 
that, whenever there is an additional interaction switched on, a typical trajectory is curved, 
its radius of curvature being large for anyons near either bosons or fermions. In this regime, 
they conclude, the mean field approximation can be justified. However, the need for an 
additional interaction to show asymmetries remains somewhat obscure. The results of Kiers 
and Weiss0 suggest that one can do without it and nevertheless find a handedness in the 
scattering: particles are defiected by a regular lattice of fiux tubes if the full coherent 
scattering is taken into account. As we shall see, these are not the only ways to make parity 
breaking manifest. 

As was already remarked in Ref. ^, broken T invariance should lead to asymmetries 
in transport coefficients. A thorough investigation of the experimentally observable con- 
sequences of T and P violation in high-T^ superconductors was carried out by Halperin, 
March-Russell and Wilczekl. The standard Onsager reciprocity relations do not apply when 
P or T are broken, and anomalous transport phenomena are expected. Working in the 
framework of an effective London Lagrangian, they find that anyon systems should exhibit 
a Hall conductivity in the absence of an applied magnetic field. 
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In this paper we set out to study a simpler model which can be expected to capture some 
essential aspects of P (and T) violation in macroscopic magnitudes. Specifically, the model 
we analyze consists of a gas of charged particles moving through a dilute random array of 
flux tubes. We shall work with a pure AB interaction, not imposing ab initio an asymmetry 
in the scattering cross-sections. We shall explicitly compute the electrical (dc) conductivity 
of the system, expecting to find a spontaneous Hall effect. From the viewpoint of mean field 
theory the graininess of the array should be unimportant. In this case one would expect to 
find a transverse conductivity similar to that of the Quantum Hall effect. 

The outline of this paper is as follows: in Section II we set up our model, in which 
flux tubes are regarded as impurities distributed randomly, each acting independently on 
electrons. After presenting the formalism for calculating transport coefficients using fleld 
theory, we use it to compute the longitudinal conductivity of the system. This is essentially 
a calculation of vertex corrections expressed in the form of a ladder equation. In Section 
HI, the parity violating, transverse conductivity is computed along these lines. Although 
divergent terms appear in the calculations, all physical quantities are flnite. Section IV 
contains a discussion of the results. 



II. LONGITUDINAL CONDUCTIVITY THROUGH AN ARRAY OF FLUX 

TUBES 

A. The Model 

We start from a Lagrangian describing a free electron gas interacting with an ensemble 
of flux tubes: 

£ = ^Fcrmi - j ■ A - p0 - pAV2m , (2.1) 

where Cpermi is the free electron Lagrangian, j and p its corresponding current and density. 

The flux tubes are located at points x^, each being the source of a vector potential {k is 
the flux per tube), 

A(x-x.) = «:^f^^^^ (2.2) 
|x - x^l^ 

and a scalar potential, 

0(x-x,) = _ (2.3) 

Then A(x) = J2xa ~ ^a), 0(x) = J2xa 'Pi^ ~ ^a)- The flux tube density will be denoted 
Hi, and the electron density Ug. 

It is known that the presence of the contact interaction, Eq. (Ol), is needed whenever 
one treats the AB problem within the framework of CS theorycHl^l. When its strength 
g is properly adjusted, the AB scattering amplitude is reproducedil. We stress that this 
interaction has nothing to do with the additional interaction introduced in Ref. ^. In our 
model, the right value for g is 

g = ±K/2m. (2.4) 
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This is half the value found in Refs. The difference is due to the fact that we are 

considering the flux tubes to be infinitely massive. Therefore the reduced mass of the 
electron/flux tube system is twice that of a pair of identical anyons. The choice of sign for 
g determines the 'hand' of the interaction. 

The vertex for interaction of flux tubes and electrons in momentum space is depicted in 
Fig. 1. The corresponding amplitude is 

K p X q 

u{p + (l,p} = g + t —. (2.5) 

m q-^ 

When g is given by Eq. (|2.4|) , this corresponds (apart from a kinematical factor) to the AB 
scattering amplitude in Born approximation. Flux tubes are assumed to be fixed objects 
with no internal excitations, so the electrons scatter from them elastically. This means that 
no frequency is carried by the interaction lines. 

As we shall show explicitly below, the 'seagull' interaction pA^ /2m can be consistently 
neglected since it contributes higher orders in perturbation theory. 

We expect that this simple model serves to approximate several many-body systems 
where some kind of AB interaction is present. Of course, more realistic interactions should 
be taken into account before it is subject to experimental verification (for a detailed study 
see Ref. ||). In this paper, it will be used as a toy model for the study of the macroscopic 
effects that we expect to arise reflecting the presence of the underlying AB interaction. 
In particular, can we make parity breaking show up anyhow in transport coefficients? Our 
answer to this question will be on the positive. 



B. Computation of the Longitudinal Conductivity 

We shall regard flux tubes as impurities randomly distributed with density << 1. 
The first effect one expects is the appearance of a finite longitudinal conductivity. This 
calculation will be done in this section, following the treatment of Ref. [l^. Having checked 
the validity of the method for this model, we shall use it in the next section to find the 
transverse conductivity. 

Consider then applying an external electric field E'^^* = —dA^^^/dt to the system. A 
longitudinal conductivity will be induced, the (linear) response being characterized by the 
conductivity tensor 

Ji = a.jEf = lua^.Af . (2.6) 

On the other hand, linear response theory yields the following relation for the current induced 
in an electron gas by an applied vector potential, 

= -f!^S,, + %(q,^))Af . (2.7) 



From ( p. 61) we see that the dc conductivity can be obtained from ( p.7|) computing to OiuS) 
with q = 0. njj(q, cu) is the (retarded) current-current correlation function. It is best 
calculated in the imaginary time formalism. This means that frequencies will be discrete 
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uJn— 27m//3 for bosons (2.8) 

= 27r(n + 1) 1(5 for fermions (2.9) 

while real frequency integrals (27r)~^ / duj are replaced by discrete sums In imaginary 

time, 

n,, (q, UJn) = f dr e^"""n,, (q, r) , (2.10) 



n,,(q,r) = -(T j|(q,r)j,(q,0)) , (2.11) 

which can be calculated using standard diagrammatic techniques. After performing the 
frequency sums one can analytically continue to real retarded time by replacing iujn 

The free electron propagator (no impurities present) is 

G%uJn,p) = {iuJn-^p)-' , (2.12) 

where = p'^/2m — ep is the electron energy above the Fermi level. 

The effect of impurities is that the electron propagator acquires an imaginary part cor- 
responding to the finite lifetime of excitations above the Fermi level. The Dyson equation 
for the propagator leads to 

GicOn: P) = . T ^— T . (2.13) 

The real part of the self-energy, Re S(a;„, p) has been absorbed in a renormalization of the 
Fermi level. 

Since flux tubes are randomly distributed, we must average over the position of each tube. 
After averaging, the Green's function in the presence of impurities becomes translationally 
invariant: (G(p, p')) = G{p)6^'^\p — p'). On the other hand, we shall be interested in 
excitations very near the Fermi surface. This will help us simplify many calculations, since 
all momenta will be strongly peaked around the Fermi value, pp- 

Our task now is to compute Im S. Diagrams with only one impurity line give (using 
translation invariance) a constant which represents a shift of the Fermi energy, irrelevant 
for our purposes. Now, under the assumptions that: a) the impurity density, n^, is low 
enough, b) scattering by impurities is weak (small k, Born approximation), then the main 
contribution to the imaginary part of the self-energy is the one shown in Fig. 2. For low 
densities the interference between scattering from different impurities is negligible, i.e., the 
scattering is incoherent. Summing over the position of the impurities gives then a factor rij. 
On the other hand, the seagull vertex, although of the same order as the diagrams in Fig. 
2, does not contribute to the imaginary part. 

These considerations lead to 



^^G'°(cUn,k)|M(k,p)|2 
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Here we have used |p|, |k| ~ pp and assumed that, for k near the Fermi surface, |M(k, p)p 
is a slowly varying function. This term is 



|n(k,p)|^ = = . (2.15) 

2m"^ 1 — cos ifi 

This is (up to a factor m/27rf^) the differential cross section da /dip for AB scatterinjll in 
the Born approximation. It is even in the scattering angle. 

Values of far from the Fermi surface give negligible contributions to the integral over 
energies. Therefore we can extend the range of integration from {—ep, +00) to (—00, +00). 
One finds 

sgn Un 

ImS(u;„,p) = , (2.16) 

Zt 



where 

= n,m I ^\u{^)\' . (2.17) 

The integrated cross section is a = [niVpT)^^. According to this, r is the mean time between 
collisions in the Born approximation (2r is the lifetime for an excitation near the Fermi 
surface). Observe that is divergent, reflecting the long range of the AB interaction. The 
integral in Eq. ( |2.17D has to be considered as regularized with a cutoff for small (f, which 



eventually has to be sent to zero. We do not write it explicitly since we shall find that it 
disappears from physical macroscopic magnitudes. 
Substitution of into ( p7[3| ) yields, 

G(^n, P) = . . (2.18) 

The computation of the polarization tensor must take into account the fact that the 
average (over positions of flux tubes) of two propagators does not equal the product of 
the separate averages. The resulting terms can be conveniently included in the form of a 
vertex term, as in Fig. 3. The vector vertex T takes account of diagrams where an impurity 
interacts with both the upper and lower electron lines in Fig. 3. Although it may seem at 
first that these terms contain higher powers in the flux tube concentration, n,, this is not 
the case since they also contain higher powers of Im S. The seagull vertex can be easily seen 
to yield a contribution of higher order in k, and therefore will be neglected. Under the same 
assumptions made above in the computation of Im E, the main contribution comes from the 
ladder diagrams (see Fig. 4), and can be found by solving the following integral equation: 

d^P' ^/ nf, / /M2 , • K qx (p-p') 



/dp ( 
— -— G{up,p')l\u{p, p')\^ + ig 



(p' X p) q X (p - p' 



|p ~ P'l^ 



(27r)2 ^ ^'^n' m |p-p'|2 

r(p',p' + q)G(a;p + a;,p' + q) . (2.19) 



(We are not writing the explicit frequency dependence of F). In general, this equation is 
very difficult to solve. However, we know that to compute the longitudinal conductivity we 
must take q = 0. The resulting equation is 
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r°(p) = 2p + I G{u„ p') \u{p, p')\'T'{p')G{u;, + p') . (2.20) 

Now, it is clear that r°(p) oc p. Write then r°(p) = (2 + A)p , where, for |p| ~ pp, A can 
be taken as independent of |p|. Eq. ( |2.20| ) gives then the following equation for A: 



A = n,m(2 + A) / ^G{u„p')G{u, + u,p') f ^\u{y,)\'^ . (2.21) 



The integral over energies can be easily performed extending the limits to (— oo, +C)o), and 
using contour integration. We quote the result since it will be used repeatedly: 

-—G[Up,p)G[Up + u,p) = ^— , (2.22) 

ZTT + i/r 

for a; > (considering ix; < does not affect our final results). Then 

A = 2^(7.)^^ + ^) > 0) , (2.23) 

with 

rf ^ = mm I ^|m((/?)|2cos¥p , (2.24) 

,;.^.--.r'^|^. (2.25) 

Ttr is the 'transport time' between collisions. The divergence in rf^ is of the same kind as the 
one in r~^, so that their difference is finite. Therefore, the dependence on the regularization 
of integrals disappears from rtr. 

Now we are ready to compute Iiij{uj) = = 0,uj). From Fig. 3(a) we read 

g2 r (Pp 



4m2/5t;'y (27r)2 
x2pi G{ujp, p)T]{p)G{up + uj,p) . (2.26) 

Eq. ( p.26| ) can now be solved. Calculations are quite standard (see Ref. [11]). The longitudinal 
conductivity is eventually found to be 

a, = ^ = 2^ . (2.27) 

Remarkably, this is a non-zero, finite quantity, in spite of the long range interaction. In 
contrast, if the interaction with impurities were Coulombian, r^^^ would be divergent, and 
therefore 0. However, the transverse nature of the AB interaction makes the divergence 

in |M((y9)p milder. 
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III. THE TRANSVERSE CONDUCTIVITY 



Here we shall use the techniques developed in the previous section to make the presence 
of a transverse conductivity manifest. 

A current perpendicular to the applied field must satisfy a relation of the form 

J, = o-xe,,Ef * . (3.1) 

Here, the presence of a transverse conductivity (J_l is a signal of P and T spontaneous 
symmetry breaking. The relation Jj = aiEf^^ already breaks T symmetry, but this is 
naturally expected since dissipative effects are present. However, T breaking is qualitatively 
different in ( |3.1| ) because the transverse current does not dissipate energy (J ■ E'^''* = 0). 

We find it convenient to choose a gauge in which the external field takes the form E*'^* = 
— V</)*^^*, or, in Fourier components, E''^* = —iq^cp^^^. Therefore 

Ji = -ta^eijqj(P'^' . (3.2) 

This means that, since Ji = — Hoi^*^^*, it will suffice to compute the density-current correla- 
tion function Hqi to 0{q), and then take the limit o; — > (Fig. 3(b)). 
The vertex r(p, p -|- q) can be expanded in powers of q as follows: 

r(p, p + q) = r°(p, p) + r (p, p + q) 

+r(-^)(p,p + q) + 0(g2) . (3.3) 

We have computed the zeroth order term in the previous section. The first order contribu- 
tion has been split into a part F*^ containing the vector q_|_ (defined to have components ^ijqj), 
and another F^'^"^-', containing q. We shall be interested only in F*^. The integral equation 
for it can be obtained from Eq. ( |2.19| ) (observe that the last term inside the brackets in Eq. 
( p.l9| ) does not contribute to F% but to F^""" ^\ since (qx p)(p' x p) = p^q-p' — (q-p)(p'-p)): 

r (p, p + q) = y ^ G{uj„ p )[-w- |p_p/|2 r (p , p ) 

+ |n(p,p')pF^(p',p' + q)|G'(a;p + ^,p') . (3.4) 



Changing the sign of g (Eq. ( p.4|) ) would reverse the flow of the transverse current. 

The ansatz that allows us to solve ( |3.4] ) is not very hard to guess. After some examination, 
one is led to write 

r^(p,p + q) = ^p(p-q±)Ai + q±A2 , (3.5) 

Pf 

where, again, Ai, A2 are independent of |p|. 

The following integrals are needed to solve Eq. (^.4[) (yj is the angle between p, p') : 



(---)p(p-q±) + -q±, (3.6) 
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— o \ ^p = — — p(p-q±) -q± • (3-7) 

Calculations are now straightforward, and lead (again for u; > 0) to 

A AA + 2)e{-ujj,)e{ujp + uj) 

' 2 UJTtr + 2 ' ^ ' ' 

A, ^ ,(A + 2)^(-.,m + .)/ _ 

Now we can use the vertex to find the part of 11 = (Iloj) that yields parity violation. 
This is (see Fig. 3(b)), 



2mP^J (27r)2 

xG{ujp, p)r(p, p + q)G{up + uj,p + q) . (3.10) 
The ansatz (|3.5| ) gives the following expression for 11'^ after integrating the angles: 

= q^l^ E( Y + A2) / |^G(o;„ p)G{u;p + u;,p). (3.11) 

Substituting the results above one finds, 

n- = iJ^^,. (3.12) 

1 — 2C(jrtr 

We have made the continuation u —itu. Highly nontrivial cancellations of t~^,Ti^ have 
concurred again to yield a finite result. Taking — > we find the transverse conductivity, 

cj^ = eV87r, (3.13) 

which is independent of Ug, Ui, and k. At first sight, this seems to lead to the nonsensical 
result that taking either or n to be zero, one still finds a finite transverse conductivity. 
Of course, this is wrong: these limits must be taken in Eq. (|3.12|) before u ^ 0, and then 
one obtains a± = 0. Also, one can check that choosing g = —K/2m, reverses the sign of the 
conductivity. 

IV. DISCUSSION 

There are two aspects of this result that deserve some explanation: first, the appearance 
of the transverse conductivity; second, its dependence on universal constants. 

It may seem somewhat surprising to find a non-zero transverse conductivity. Consider 
the motion of a single electron in the presence of a flux tube. The differential scattering 
cross section ( ^.15] ) is invariant under ip — > —ip, i.e., a typical electron trajectory is not 



deflected. In order to make asymmetries appear, March-Russell and WilczekO were forced 
to introduce an additional interaction (in the form of a generic phase shift). Then parity is 
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broken through interference terms. This is exphcitly illustrated by Suzuki et al. in Ref. [T^ , 
where the additional interaction is taken to be a hard-disk repulsion. The situation here is 
different, since the interaction with the flux tubes is purely AB. Of course, one could argue 
that we actually have another interaction, namely that with the external electromagnetic 
field. But we have reasons to believe that this by itself is not the origin of parity breaking. 
We have analyzed a closely related problem: the bremsstrahlung for emission of soft photons 
in the presence of a flux tube. Vertex corrections have been taken into account. It turns 
out that, after summing over the scattering angle of the emitted photon, the corresponding 
cross section for electrons is invariant under parity. It seems that an electron interacting 
with a flux tube emits photons in a left-right symmetric way. 

There is however another difference between our problem and that of the scattering of a 
single electron from a single flux tube: it is essential in our calculation to take into account 
the many-body effects. Appart from yielding a damping of excited electrons, many-body 
effects are present in the ladder resummation that we perform to find the vertices. One can 
check that if only the first term in the series is included, the transverse conductivity turns 
out to be zero. Of course, it is not consistent to do so, as we have argued before. We need 
to sum all the terms and it is then that the asymmetry appears. Therefore parity breaking 
does arise at the macroscopic level, although its microscopic origin is hidden. 

The universal dependence of the transverse conductivity is less surprising once it is 
noticed that it is of the kind familiar from the Quantum Hall effect. This seems to provide 
support to the mean field approximation, since we find a behaviour analogous to that of an 
electron gas in the presence of a magnetic field. However, we are dealing with a perturbative 
approximation, so we do not expect to find the Landau level structure responsible of the 
integer-spacing in the Hall conductivity. Nor, evidently, can our model take into account 
the null longitudinal resistivity that appears simultaneously with it. Nevertheless, it is 
remarkable to find in this simple model such a semi- quantitative agreement with mean field 
theory. 

The procedure we have described above can be modified to compute other transport 
coefficients, such as thermal conductivities; these will also show parity violation. However, 
it seems more interesting to take one step further and consider the conductivity of a gas of 
anyons with the full CS interaction taken into account. Calculations with this model are 
much more difficult, but the outcome would be a full description of transport phenomena in 
anyonic systems. 
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FIGURES 



FIG. 1. Vertex for interaction between an electron (continuous line) and a flux tube (square). 
The amplitude is n(p + q, p) = 5 + i^^^- 

FIG. 2. The main contribution to the imaginary part of the self-energy. The impurity lines 
give a contribution |u(k, p)p. 

FIG. 3. (a) The current-current correlation function Hij. (b) The parity violating den- 
sity-current correlation function IIq^. The electron lines correspond to the propagator G(u;„,p) 
dressed by the interaction with impurities (fig. 2). 



FIG. 4. Ladder equation for the vector vertex T. 
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